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ABSTRACT 

A number of recent lines of evidence point towards the presence of hot, out- 
flowing plasma from the central regions of compact Galactic and extragalactic 
X-ray sources. Additionally, it has long been noted that many of these sources 
exhibit an “excess” continuum component, above ~ 10 keV, usually attribut- 
ed to Compton Reflection from a static medium. Motivated by these facts, as 
well as by recent observational constraints on the Compton reflection models — 
specifically apparently discrepant variability timescales for line and continuum 
components in some cases - we consider possible of effects of out-flowing plasma 
on the high-energy continuum spectra of accretion powered compact objects. We 
present a general formulation for photon downscattering diffusion which includes 
recoil and Comptonization effects due to divergence of the flow. We then develop 
an analytical theory for the spectral formation in such systems that allows us 
to derive formulae for the emergent spectrum. Finally we perform the analyt- 
ical model fitting on several Galactic X-ray binaries. Objects which have been 
modeled with high-covering-fraction Compton reflectors, such as GS1353-64 are 
included in our analysis. In addition, Cyg X-3, is which is widely believed to 
be characterized by dense circumstellax winds with temperature of order 10 6 K, 
provides an interesting test case. Data from INTEGRAL and RXTE covering the 
~ 3 — 300 keV range are used in our analysis. We further consider the possibility 
that the widely noted distortion of the power-law continuum above 10 keV may 
in some cases be explained by these spectral softening effects. 
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1. INTRODUCTION 

Recent observational and theoretical evidence suggest that accretion-powered X-ray- 
sources, of both the Galactic and extragalactic variety, may exhibit outflowing plasma, i.e. 
winds, emanating from a compact- region near the central source (e.g. Elvis 2003; Arav 2003; 
Brandt & Schulz 2000; Proga & Kallman 2002). Comptonization effects in those putative 
outflows are likely to alter the intrinsic continuum spectra of accretion powered compact 
objects. The basic idea is that electron scattering of photons from a central source entering 
the expanding outflow experience a decrease in energy (downscattering). The magnitude of 
this decrease is of first order in v/c and in E/m e c 2 where v is the outflow speed, c is the 
speed of light, E is the initial photon energy, and m e is the electron rest mass. 

The basic idea is depicted in Figure 1. There present a simple explanation of the 
diverging flow effect on the photon propagation through the medium. A photon emitted 
outwards near inner boundary and then scattered at a certain point by an electron moving 
with velocity v lt is received by an electron moving with velocity V 2 as shown with frequency 
i/ 2 = [1 -p (vj — V 2 ) • n/c] where n is a unit vector along the path of the photon at the 

scattering point. In a diverging flow (vj — V 2 ) • n/c < 0 and photons are successively red- 
shifted, until scattered to an observer at infinity. The color of photon path (in Figure 1) 
indicates the frequency shift in the rest frame of the receiver (electron or the Earth observ- 
er). On the other hand, referring to the right-hand side of Figure 1, in a converging flow 
(vj — v 2 ) • n/c > 0 and photons axe blueshifted. 

The classical Compton (recoil) effect < A E > jE ~ —E/m e c 2 has been well understood 
for a long time. Basko, Sunyaev h Titarchuk (1974), hereafter BST74 first studied the down- 
scattering effects in the interaction of X-ray radiation of the central source with the relatively 
cold atmosphere of the optical companion in a binary system. They predicted the shape of 
the X-ray reflection spectrum of the companion and applied these results to the Her X-l 
system. X-ray observations of Her X-l by a number of groups (e.g. Sheffer et al. 1992, Still 
et al. 2000 and etc.) confirmed this prediction. In particular, the numerical calculations by 
BST74 demonstrated that the reflection of the bremsstrahlung spectrum distorted the high- 
energy continuum above 10 keV producing the characteristic continuum excess feature, or 
“bump” , in the > 10 keV spectrum. This feature was also identified by Sunyaev & Titarchuk 
(1980), hereafter ST80, who found that the transmitted, downscattered spectrum is formed 
as a result of the reprocessing of X-ray radiation from the central source in an ambient 
spherical cloud. From these facts one can conclude that this distortion of the continuum 
is not an intrinsic property of the particular incident photon distribution, but it is rather 
a result of down-scattering effects from X-ray photons propagating through the intervening 
medium. 
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In this paper we show that this is precisely the case. We formulate the problem of the 
photon diffusion in generic terms and demonstrate a specific solution for bulk motion where 
effects of a recoil and divergence of the flow are taken into account. The emergent spectrum 
as an outcome of this solution has all these aforementioned features of the downscattering. 

The details of the radiative transfer problem taking into account down-scattering effects 
(recoil and Doppler effect in the divergent flow) and its solution are given in §2. In section 
3, we apply our model to observational data for several Galactic X-ray binaries. We attempt 
to demonstrate the Compton down-scattering effects in relatively cold outflows for which 
temperature is of order 10 2 * * * 6 K and disk atmospheres which are necessary constituents of X- 
ray sources. For example, Cygnus X-3 is believed to be characterized by dense circum-stellar 
(or circum-disk) winds, provides an interesting test case. A number of additional sources, 
which have been noted in the literature to exhibit evidence for strong “Compton reflection”, 
e.g. GS1353-64 and GX339-4, have also been included in our analysis. We consider the 
possibility that this well documented distortion of the power-law continuum above 10 keV, 
may in some cases be due to downscattering from out-flowing plasma rather than from a 
static reflecting media. We further speculate that this could possibly apply to extra-galactic 
(i.e. AGN) as well as Galactic accretion objects (§4), although we do not present any analysis 
of AGN here. A summary and conclusions follow in §5. 


2. RADIATIVE TRANSFER IN A BULK OUTFLOW 

2.1. The Basic Downscattering Problem 

The problem of photon propagation in a fluid in bulk motion has been studied in detail in 
a number of papers [see e.g. Blandford & Payne (1981), hereafter BP81; Payne &: Blandford 

(1981), hereafter PB81; Nobili, Turolla & Zampieri (1993); Titarchuk, Mastichiadis &: Kylafis 

(1997); Titarchuk, Kazanas & Becker (2003, hereafter TKB03); Laurent & Titarchuk (1999), 
(2001), (2004, hereafter LT04]. In particular, TKB03 present a general formulation and a 
solution of the spectral formation in the diverging outflow. They demonstrated that the 
resulting spectrum can be formed as a convolution of energy and space diffusion solutions. 
The spread function, given by a Green’s function formulation, applied to a monochromatic 
injection yields a redshift-skewed line which is a power law at low energies. Furthermore, 
Monte-Carlo simulations show that the line width has a strong dependence on u/c, optical 
depth r and the energy of the monochromatic line Eq only when the flow temperature kT e < 1 
keV (LT04). TKB03 and LT04 also establish that the shape of the spectral line should be 
closely related to the photon source distribution in the flow. Comparison of the theoretical 
spectra with data show that observed red-skewed lines are formed as a result of transmission 
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of the X-ray radiation through outflows of moderate Thomson optical depth (r 0 > 1). 

To apply these Radiative Transfer results to observational data, we have developed a 
generic analytical formulation which leads to a simple analytic expression for modification 
of the emergent spectrum due to recoil and velocity divergence effects in the flow. For the 
recoil effect, we extend the results of ST80 who show that the down-scattering feature, or 
’’bump”, can appear superposed on power-law spectra with indices a < 1 as a result of 
photon diffusion through a static cloud. From ST80 and calculations presented here, one 
might conclude such bumps are not a necessarily a feature of disk reflection, but could 
instead be a generic feature of the photon reprocessing in a relatively cool, ambient plasma 
characterized by temperatures of order 10 6 K. 


2.2. Main Equations and Solution of Downscattering Problem 

Let N(r ) = iVo( r o/' r ) b be the radial number density profile of an outflow and let its 
radial outward speed be 


v b /c= (M ou t/47rcN 0 rl)(r 0 /r) 2 b = u 0 (r 0 /r) 2 6 (1) 

obtained from mass conservation in a spherical geometry (here M ou t = 47r r 2 v b N). The 
Thomson optical depth of the flow from some radius r to infinity is given by 


CO 

N e (r)a r dr = arNor^ro/r)^ 1 j (b - 1), (2) 

where N e (r) — N(r ) is the electron density, cr T is the Thomson cross section, r 0 is a radius at 
the base of the outflow. Because our final results are independent of the velocity and density 
profiles below for simplicity of presentation we use b = 2 (a constant velocity outflow). In 
this case To = r(r 0 ) = a^Noro. 

The transfer of radiation within the flow in space and energy is governed by the photon 
kinetic equation (BP81, Eq. 18) for the photon occupation number n(r,i/), which in steady 
state reads 



v b /c 11 1 V • (v b /c) dn 1 d(z 4 n) 

' Vn+ s v fr Vn)+ - 3 —^— z Tz + 7~ 




( 3 ) 


where z = hu/m e c 2 is the dimensionless photon energy, k = N e (r)a T is the inverse of the 
scattering mean free path, v b = v b e r , is the flow velocity, e r is the radial unit vector and 
j{r, v) is the photon source term. 
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The spectral flux F(r, i /) (PB81) is given in terms of n(r, u) by 

. 1 _ 1 dn 

FM= ~W) 3 ViV ^' () 

and must satisfy the following boundary conditions: (a) Conservation of the frequency inte- 
grated flux over the outer boundary, namely 

J F(r, v)dv oc r~ 2 as r — + oo. (5) 

(b) The occupation number should be zero at the inner boundary r = r 0 , where the bulk 
outflow starts and the outflow radiation is zero, i.e. 


n(r 0 ) = 0. (6) 

It is worth noting that in TKB03 the goal was to demonstrate the pure effect of the flow 
divergence on a spectral line as a result of multiple scatterings with the plasma free electrons. 
They thus omitted the recoil term in the left hand side of equation (3). They were also 
interested to applying their solution to the spectral formation of iron lines for which photon 
energies are about 6.4 keV and lower. For such energies the recoil effect can be safely 
neglected (see ST80). On the other hand if one is interested in the modification of continuum 
spectrum at energies higher than 10 keV by multiple down-scattering events in the “cold” 
medium this effect must be taken into account. Below we show how one can generalize the 
TKB03 solution including the recoil term in the radiative transfer equation. We apply a 
particular method for the separation of variables suggested in Titarchuk (1994), hereafter 
T94, and successfully used in TKB03, to obtain our solution. 

In order to further proceed with the solution derivation we should note that Laurent 
and Titarchuk (2004) find that the coefficient in the diverging term of equation (3) e = 
V • (V b /c)//c can be replaced by a constant i. Strictly speaking e is a function of r and 
consequently of r, namely e = 2(u&/c)/r. 

Using the Fokker-Planck equation [see Eq. (3) in TKB03] and a method developed by 
Titarchuk, Masticbiadis & Kylafis (1997) (see Appendix D there) one can relate the mean 
energy change per scattering, < A E >, for photons undergoing numerous scatterings in the 
flow to /? = v/c: 

<AE>zt-[fVv/{cK)]E^-(2f(3/T 0 )E (7) 

where v = ve r is the flow velocity, e r is the radial unit vector, k = N e a-r is the inverse of 
scattering mean-free-path l. The numerical factor /, in formula (7) is of order unity and 
LT04 obtain its precise value using their MC simulations. 
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In LT04 (specifically, Figure 2 of that paper) the escape photon distribution ip(t) for 5 
energy bands is present. In the plot the time is given in light crossing time units t cross = A r/c 
where A r is the outflow cloud thickness. Photons which escape without any scattering are 
at 6.6 keV. The model parameters kT e = 0.1 keV, r 0 = 4, /? = v/c = 0.1. 

The simulated time distribution is fitted by an exponential <p(t) = Cm exp(-at/t cross ), 
with a equal to 0.67 identical to that predicted by diffusion theory [see Sunyaev & Titarchuk 
in 1985 (hereafter ST85)]. ST85 show that the average number of scatterings N av in the 
shell of optical depth r 0 is 3 tq / 8 so the average photon scattering time is t av = N av l/c = 
3 to t cross /8, where l = 1 /Nq<j t = A r/r 0 . ST85 also show that the time distribution for 
scattered photons in any bounded medium is an exponential, ip(t) = Cm exp(—t/t av ) where 
C N = l/t av . 1 ). 

In this particular case /(f) = Cm exp(— t/t av ) = Cm exp[— (2/3)/(ro/4)(f/t cross )], pre- 
cisely what is obtained in the LT04 simulations. Thus, the average energy of photons escap- 
ing after N av -scatterings is < E > sc ~ (1+ < A E >) Nav . Using formula (7) one can obtain 
that 

< E > sc « (1 - 2fP/r 0 ) N "-E 0 . ( 8 ) 

For the particular case of r 0 = 4 and (3 = 0.1 LT04 find that the original photon energy 
Eq = 6.6 keV is reduced to Eq = 5.4 keV, after N av —scatterings in the outflow [the emergent 
spectrum is shown in Fig 1, (LT04) left panel]. The analytical estimate obtained from 
formula (8) is close to this value of < E > sc = 5.4 keV for / ss 2/3. Thus one considers an 
approximation in which the diverging term in equation (3) is independent of space variable. 


2.3. Downscattering solution: Emergent spectrum 

According to a theorem (T94, appendix A), the solution of any equation whose LHS 
operator acting on the unknown function n(r, v) is the sum of two operators L r and L u , which 
depend correspondingly only on space and energy and the RHS, j(r,u), can be factorized, 
i.e. 

L r n + L^n = —j (r, v) = -f{r)<p( u). (9) 

The boundary conditions independent of the energy is, i.e. 

L^n = 0 as r-+oo, L^n = 0 for r = r 0 , (10) 

n is given by the convolution of the solutions of the time-dependent problem of each operator, 
namely 

roc 

n(r,u) = / P(r,u)X(v,u)du. 

Jo 


( 11 ) 
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Above w is the dimensionless time (which can be the Thomson dimesionless time «t = iV 0 cr x ct 
depending on the specific forms of operators L r and L v ) and P(r, it) is the solution of the 
initial value problem of the spatial operator L r 

f£ = irP, P(r,0) = /(r) (12) 

with boundary conditions 

L^P = 0 as r — ► oo, L^P = 0 at r = ro, (13) 


and X(u, u) the solution of the initial value problem of the energy operator L u 



^ = L„X, X(z,0) = tp(z) 

(14) 

with boundary conditions 


z 3 X — + 0 when z —* 0, oo. 

(15) 

Thus we have 

L r P = Vb / C - VP + 1 V • {-VP). 
K 3 KK 

(16) 


Following, the LT04 arguments we replace the coefficient of the diverging term by a 
constant e = 2q(vb/c)/TTfl 1 where a numerical factor q is order of unity (see details in 
the end of §2.2). 

In the case where the photon energy is due to a recoil effect and Comptonization effects 
in the diverging flow, we have 


dX T v 1 8X 1 d(z 4 X ) 

aT = L - x = 3“ & + 

x(z, 0) = ip(z)/z 3 


with boundary conditions 


(17) 

(18) 


z z X — * 0 when z — ► 0, oo. (19) 

We transform Eq.(17) introducing a new unknown function Y = e^ 4c ^ u z 4 X for which equa- 
tion has a form 



( 20 ) 


The problem for equation (20) with appropriate initial condition Y(z,0) = zip(z) and 
boundary conditions Y — > 0 when z — ► 0, oo is an initial value problem for the first order 
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partial differential equation and it can be found using the method of characteristics (see e.g. 
TKB03). The differential equation for the characteristics is 

dz 

= 

that solution is 

u — -[ln(z + e/3 )jz - ln(z 0 + e/3)/z 0 )], 

where zo is the dimensionless energy at u = 0. Because Y (zo, 0) = zo<p(zo) 
the characteristics the solution of the problem (17-19) is 

J(z,u ) = z z X(z,u) = e~ ieu/z ip[ip e (z,u)]/[z'tp e {z,u)}, (23) 

where z 0 = if e {z, u ) = (<e/3)/[(l +e/3 z) exp(— eu/3) - 1] is found from Eq. (22). Substitution 
of J(z,u ) from Eq.(23) into Eq.(ll) gives us the emergent spectral shape 

1 rV'max,e{z) 

F E (z,e) = [t~ 2 F(t,z)} | r ^ 0 oc- / e _4e “ /3 [^ e (z,u)]~V[^e(A«)]^(u)d« (24) 

Z Jo 

where u max>e {z ) = (3/e)ln(l + e/3 z) > 1 and V{u) oc [t _2 <9P/<9t[t,u] | t ^ 0 using the ex- 
pressions for F(t,z ) and P(r, u) [see Eq. (4), Eqs. (12-13) for the definition of P(r, z) and 
P(t,u) respectively]. 

This result is a generalization of ST80 and TKB03 results. ST80 derived the spectra 
when the down-scattering effects due to the recoil was taken into account. On the other 
hand TKB03 considered the downscattered spectra as a result of the divergence of the flow. 
For example, one can obtain the ST80 formula (36) for the recoil spectrum assuming e — > 0 
in Eq. (24). In fact, 

lim ip(z,u) = (1/z-u)- 1 

e—*0 

and u maXt o(z) = 1/z. 

That formula (24) for e = 0 ( the recoil case ) is generic and it is valid for any geometric 
configuration of the plasma cloud, e.g. a disk, or a spherical cloud, as well as for any photon 
source distribution within the cloud, e.g. uniform, or central illumination distributions. In 
fact, what we use here to derive this formula are the particular properties of the diffusion 
operator of the left hand side of the main equation (9), the boundary conditions (10) and 
the source term in right hand side of equation (9). Namely i. the diffusion operator should 
be a sum of two operators, one is the space diffusion operator and another one is the energy 
operator, ii. the boundary condition are independent of energy and iii. the source function 
is factorized [or presented as a linear superposition of the products of f(r)p(v)}. All these 
mathematical properties are generic for the diffusion problem and independent of any specific 


( 22 ) 

is conserved along 
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geometry and spectral and space source distribution (see ST80 and T94 for more details and 
particular examples). 

Below we demonstrate how formula (24) can be simplified by exploiting the fact that 
downscattering energy change AE/E is proportional to e and E/m e c 2 . 

It should be noted that the down-scattering modification of the spectrum occurs when 
photons undergo multiple scatterings u. It is always on the order of the average number 
of scatterings N av — J 0 °° uP(u)du/ J 0 °° V(u)du. Thus we can expand the integrand function 
W e (z,u) — e~ 4EU ^[ij) £ (z,u)]~ 1 if[tj} e (z,u)] over u in formula (24): 

W(z, u) e w W E (z, 0) + W'(z, 0 )u. (25) 

This expansion is valid because we consider the case when N av -C 1/s (r Tt0 ll, £ < 1). 
Substitution of Eq (25) into Eq. (24) gives us 

f(v) oc z- l [W e {z, 0) + N av W' e {z, 0)], (26) 

where 0) is the incident spectnnn <p(z) and 

W e '(z,0)]/z = <p(z) |{ln[z 0 y>(zo)]} / zo | u= o^iu=o - 4e/3 j , (27) 

where z 0 = ^ E (z,u) = (e/3)/[(l + e/3z)exp(— eu/3) — 1] and 

£i— < 28 > 

Substitution of equations (27, 28) and <^{z) — z~ l W e (z 1 0) into (26) leads us to the formula: 

E e (z, e) oc c^(z) { 1 - (4e/3 )N av + iV ai) z 2 (l + e/3z)[hi(zo<p(zo))](J„=o}- (29) 

The most interesting case is the down-scattering modification of Comptonization spectrum 
(see e.g. ST80 and Titarchuk 1994) that can be well fitted by a powerlaw spectrum with an 
exponential cutoff, namely 

<Pcomp(z) ~ exp (— z/z*). (30) 

The cutoff energy E, is related to the Compton cloud electron temperature kT e , namely 
E, se 2 kT e . For <p comp (z), it is evident that 

[ln(z 0 ^(z 0 ))y u= o = (1 - a)/z ~ l/z* (31) 

Using equation (31) we transform formula (29) as follows 

Ee{z, e) oc ip comp {z){ 1 - (4c/3)7V ot . + N av z( 1 + e/3z)[(l - a) - zjz *]}. (32) 
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The modification of the absolute normalization of the spectrum due to the down-scattering 
effects is quite obvious — the relative change of the normalization is 1 — (1 + a)eN av / 3 . 
Thus we can rewrite Eq. (32) as the final formula of the spectral shape as follows 

Fe(z,£) cx <Pcomp{z){ 1 + N av z[( 1 - a) - e/3z* - z/z t ], (33) 

where N av = N av /[ 1 — (1 + a)eN a v/ 3]- The second term in parenthesis of formula (33) 
describes the pile up and softening of the Comptonization spectrum <p C om P {z) due to the 
down-scattering effect in the outflow. In Figure 2a we present ETe diagram for various N av 
and e. The down-scattering bump and softening of spectrum are clearly seen in this plot 
[see also Fig. 2 where a ratio of the Comptonization models to the incident spectrum (30) 
is plotted]. 

Below we apply formula (33) to the X-ray spectral data for a number of sources. The 
self-consistency of application of this formula for data fitting can be checked by comparison 
of the best-fit parameter N av and u mox e (z max ) [the formula for u maxe (z) is shown just after 
Eq. (24)]. Our inferred best-fit parameter N av -C u maXt g(z Tnax ) for all fits and thus our 
formula (33) as the analytical approximation of (24) is valid for all cases considered. 


3. DATA ANALYSIS 

To test these ideas we obtained public data from the INTEGRAL and HEASARC 
(RXTE) archives for a sample of Galactic X-ray binaries (Table 1). The Cygnus region 
was covered extensively by INTEGRAL in the early stages of that mission, Performance 
Verification Phase, (November-December 2002), subsets of which are in the public domain 
at the time of this submission. Specifically, we selected coverage of Cyg X-3 and Cyg X-l 
obtained during revolutions 22-25. Statistics are limited by the small amount of useful data, 
given the early problems with data gaps. We extracted spectra from JEM-X and SPI, using 
the OSA release version 3 software. Subsequent model fitting was done with the XSPEC 
spectral analysis package (version 12.10a:) modified to include a local implementation of the 
model described in section 2. We also analyzed the JEM-X and SPI data for Cyg X-l, 
applying the same model. For comparison, we also obtained and analyzed contemporaneous 
RXTE data for both sources. In all cases, we used the RXTE calibration database and 
software release current as of April, 2004. We also utilized RXTE data for the additional 
sources GX 339-4 and GS 1354-63, each of which have been modeled by various other groups 
- e.g. Gilfanov et al (1999); Nowak, Wilms & Dove (2002); Pottschmidt et al (2003); Vilhu 
et al (2003) - as a power-law- exponential plus a Compton ’’reflection” continuum (see e.g. 
Magdziarz & Zdziarski 1995). For example, GS 1353-64 was found to require a large reflection 
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continuum, R ~ 0.3 — 0.6 (Gilfanov et al 2003). In addition to the Comptonized continuum, 
iron line structure - a gaussian line component and absorption edge feature - was found 
to improve the quality of fit at low energies. In addition, a systematic error component 
of 1% was added to the PCA data. We note that there is a significant cross-calibration 
discrepancy between the INTEGRAL SPI and JEM-X instruments (e.g. Paizis et al 2003). 
We have assumed that the absolute flux calibration of SPI is reliable (Attie’ et al. 2003; 
Sturner et al. 2003), and renormalized the JEM-X model fits accordingly. 

In Table 1 we list our source sample, and the inferred parameters of the out-flow- 
Comptonization from our analysis. Here T = a + 1 is the (photon) spectral index. N av is 
the average number of scatterings experienced by a typical photon in the outflow (section 
2). The origin of the data RXTE (PCA plus HEXTE) or INTEGRAL (JEM-X plus SPI) 
are also indicated. We note that two Cyg X-3 observations represent very different intensity 
states for that source. For that matter, the GX 339-4 and GS 1353-64 observations used 
correspond to high-intensity states for each of those sources (but in case, the spectral energy 
distributions represent the low-hard state). 

Figure 3 further illustrates the down-scattering effect. Plotted there is our best fit model 
curve, using GX339-4 as a test case, compared to a simple exponentially-folded power-law 
fitted to the same data, with an accompanying residual plot. T his illustrates the improvement 
in the fit above 10 keV. It is clear from these plot that the down-scattering effects lead to an 
improved fit above ~ 10 keV. Also one can see the photoabsorption features below 10 keV 
and higher K-edge (in 7-9 keV energy band) in the residual plot (Fig. 3a). It is worth the 
similar edge features along with the strong K a lines are detected during X-ray superbursts 
(see Figs. 5 and 9 in Strohmayer & Brown 2002). There is a high probability that they are 
originated in the radiation driven outflows during burst events. 

In Figure 4, we show an example of one of our model fits, in this case to GS 1353-64 
(Figure 5 is the same result, but plotted in photon space). 


4. DISCUSSION 

We have shown that down-scattering modification of the primary photon spectrum by 
an outflowing plasma is a possible mechanism for producing the continuum excess in the 
~ 10 keV spectral region. This is usually attributed to Comptonization by a static reflector, 
such as a downward- or obliquely-illuminated accretion disk, although the overall continuum 
form differs from that of the basic Compton reflection form. We thus suggest, that in at least 
some cases, the outflow down-scattering effect rather than the standard Compton reflection 
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mechanism is responsible for the observed ” excess “ hard- X-ray continuum. 

It is reasonable to expect that outflowing plasma, in the form of stellar or putative 
disk winds will effect the emergent spectra of compact binaries. Collimated outflows are 
well known in certain objects, the so called Galactic micro-quasars, and there is strong un- 
collimated or weakly- collimated outflow in additional objects; as noted Cyg X-3, where in 
fact spatially extended emission has been resolved (Heindl et al. 2003), and of course the 
stellar winds in Cyg X-l have been extensively studied. Disk winds may have been observed 
directly in Circinus X-l, e.g. Brandt & Schulz (2000). Recently, evidence for disk winds 
in AGN has emerged (e.g. Elvis 2001; Arav 2003). Further observational confirmation is 
needed, but if present, these winds could produce the down-scattering effects we propose. 

In addition to evidence for outflows, recent observations suggest that the putative reflec- 
tor is, in at least certain objects, necessarily much larger than any reasonable disk size(e.g. 
Mattson & Weaver 2004; Markowitz, Edelson & Vaughan, 2003). These arguments are based 
on temporal signatures; specifically the lack of a prompt response of the ” reflected” emis- 
sion to the continuum assuming light travel times within reasonable accretion disk spatial 
scales. While an accretion disk torus has been suggested as an alternative, i.e. more remote, 
reflector, the same temporal signatures could be reproduced within the context of a central 
source and an ambient outflow. 

We note the concerns of some authors, e.g. Miller et al. (2004), that the hot disk 
inner-disk component and high frequency QPOs ( few times 100 Hz) seen in accreting Black 
Hole (BH) at the highest inferred accretion rates would not be visible through an “optically 
thick” outflow. Those authors thus rule out any possibility of ambient spectral reprocessing 
and distortion of the iron line feature, such as broadening and red-wing enhancement, by 
the outflow. On the other hand, Laurent & Titarchuk (2004) infer the outflow Thomson 
optical depth r 0 of the out-flowing medium from parameter fitting using XMM and ASCA 
measurements of broad, red-shifted iron lines. They find that r 0 never exceeds 2-3 in 
any of the cases analyzed. Presumably, the observer sees the radiation of the BH central 
source through the ”haze” of the moderate optical depth. Furthermore Titarchuk, Cui & 
Wood (2002), hereafter TCW02 give a precise model for the loss of the modulation due 
photon scattering. It follows from TGW02 that the outflow optical depth r ro would need 
to be around 16 and higher in order to suppress QPO amplitude of frequency 100 Hz [see 
formula (5) in TCW02], We have shown that this very optically thick outflow is ruled out by 
observations. Our results confirm the LT04 results in the sense that all effects of reprocessing 
and line distortion can plausibly occur in outflows characterized by moderate optical depth. 

Furthermore, recently Laming & Titarchuk (2004), hereafter LaT04, formulate a gener- 
ic problem of the outflow illumination by the hard radiation of the central object and they 
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calculate the outflow temperature and ionization balance as a function of the ionization pa- 
rameter. Natural assumptions regarding the X-ray spectral distribution of the central source 
radiation (Comptonization like spectra) and velocity distribution in the outflow (constant 
velocity wind) are applied. They find that iron K a photons are generated by the absorption 
of X-ray photons at energies higher than the K-edge (i.e > 7 keV) (see these K-edge features 
in Fig. 3). Electron scattering of the K a photons within the highly ionized expanding flow 
leads to a decrease of their energy (redshift) which is of first order in vfc (this is clearly 
illustrated in Fig.l) This photon redshift is an intrinsic property of any outflow for which 
divergence is positive. LaT04 find the range of the parameter Inner Radius/Z/40 (which is 
proportional to the inverse of the so called “ionization parameter” used in the literature) is 
about 10 13 cm and the range of the wind temperature is about 10 6 K when the observed K a 
lines are produced in the wind (where L40 is the source luminosity in 10 40 ergs s -1 ). They 
also find that the equivalent widths of red skewed Fe Ka originated in the wind is order of 
keV. 


5. CONCLUSIONS 

We have developed an analytic formulation for the emergent spectrum resulting from 
photon diffusion in a spherically expanding Comptonizing media, characterized by two model 
parameters: an average number of scatterings in the medium N av and efficiency the energy 
loss in the divergent flow e. Application of the model to high-energy spectra of several com- 
pact binaries, which have been previously modeled within the Compton reflection scenarios, 
seems to lead to a satisfactory representation of the data. We suggest that in some instances, 
the apparent ~ 10— keV continuum enhancement seen in Galactic and extra- galactic sources 
is due to down-scattering effects associated with an outflow rather than to reflection by a 
disk. 


We may also conclude that scattering and absorption of the primary line photons in the 
relatively “cold” outflow (that temperature is few times 10 6 K) lead to the down- scattering 
modification of the continuum and to the formation of red- skewed lines (that is a more natural 
and probable mechanism than the general relativistic effects in the innermost part of the 
accretion flow). 

In future work, we will employ detailed Monte-Carlo calculations to further explore the 
range of validity. In addition to Galactic binaries, extra-galactic sources may be characterized 
by out-flowing plasma. With the expanding database of INTEGRAL observations, which 
will hopefully lead to an increasingly accurate characterization of the high-energy continua, 
it may be possible to to explore this idea. In addition to the continuum, we will explore the 
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possibility of the line feature formation in the wind. 
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Table 1. Source Sample and Comptonization Model Parameters 


Source ID 

Observatory 

r 

N av 


6 

Flux 

xl 

DoF 

Cyg X-l 

INTEGRAL 

1.74 

0.8 

175 

0.12 

10.01 

1.544 

1359 

Cyg X-l 

RXTE 

1.51 

1.7 

219 

0.12 

17.1 

1.651 

80 

GX339-4 

RXTE 

1.63 

3.2 

299 

0.10 

4.27 

0.988 

123 

GS1353-64 

RXTE 

1.35 

1.9 

103 

0.10 

5.2 

1.015 

122 

Cyg X-3 

INTEGRAL 

2.61 

1.6 

160 

0.10 

5.57 

1.218 

1550 

Cyg X-3 

RXTE 

2.03 

1.01 

270 

0.10 

13.1 

1.068 

68 


Flux is for 3-100 keV, units of 10 9 ergs cm 2 s 1 
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XT 

Fig. 1. — On the left side: Schematic diagram depicting the wind geometry. The outflow 
(wind) originates at the inner radius. The electron optical depth of the wind is of order 
unity. A photon emitted near the inner boundary and subsequently scattered by an electron 
moving with velocity v x , impinges on an electron moving with velocity v 2 as shown. The 
change in frequency is i/ 2 = tq [1 + (vj — v 2 ] • n/c) where n is a unit vector along the path of 
the photon at the scattering point. In a diverging flow (v! — v 2 ) • n/c < 0 and photons are 
successively redshifted, until scattered to an observer at infinity. The color of photon path 
indicates the frequency shift in the rest frame of the receiver (electron or the Earth observer). 
On the right side: In a converging flow (vj — v 2 ) • n/c > 0 and photons are blueshifted. 
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Fig. 2. — a) ETe diagram representing the thermal Comptonization spectra that result 
from photon propagation through a divergent outflow. The photon source is in the center 
of the outflow shell. Solid line: the emission of the central source for which a spectrum is 
E~ a exp(—E/E») where a = 0.5 and F* = 50 keV. Dashed, dot-dashed and dotted lines: the 
emission not escaping radiation from the outflow shell for that average number of scattering 
N av = 2, 3, 5 respectively. The effective coefficient of the outflow divergence e = 0.2. b) 
Similar to Figure 2a, except that a ratio of the Comptonization models to the exponentially 
cutoff power-law spectrum is plotted. 
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Fig. 3. — a) Comparison between our Comptonizaton model (solid curve), and a simple cutoff 
power-law form, i.e. ~ KE~ r exp(—E/ E q) (dashed curve). Both curves resulted from fits to 
the overlayed data points, from GX 339-4. No additional parameters (such as absorption or 
emission due to Fe) were included. It is evident that the Comptonization model represents 
the data more accurately above ~ 10 keV (the improvement in xl is about 10%). As detailed 
in the text, we interpret the ’’excess” flux in that regime as being due to down-scattering of 
hard photons in a diverging outflow, b) Fit residuals (for GX 339-4) above 10 keV for the 
cut-off powerlaw model (triangles) and our out-flow Comptonization model (filled circles). 
As evidenced from this figure, the down-scattering effects provide a viable explanation of the 
10-100 keV ’’excess” continuum. 
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Fig. 4. — Example of the the application of our model to observational data. In this case, 
flaring state RXTE data for GS 1351-64 are shown. This source has been noted by others 
as a strong “Compton reflection” source among the Galactic X-ray Binaries. The count rate 
data, the folded model, and the residual are plotted. In this case, the chi-square per degree 
of freedom was xl = 102. 




